Generation of electromagnetic fields in string cosmology 
with a massive scalar field on the anti D-brane 
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We study the generation of electromagnetic fields in a string-inspired scenario associated with a 
rolling massive scalar field <j> on the anti-D3 branes of KKLT de Sitter vacua. The 4-dimensional DBI 
type effective action naturally gives rise to the coupling between the gauge fields and the inflaton <j>, 
which leads to the production of cosmological magnetic fields during inflation due to the breaking 
of conformal invariance. We find that the amplitude of magnetic fields at decoupling epoch can be 
larger than the limiting seed value required for the galactic dynamo. We also discuss the mechanism 
of reheating in our scenario and show that gauge fields are sufficiently enhanced for the modes deep 
inside the Hubble radius with an energy density greater than that of the inflaton. 
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Introduction - The development of string theory has 
continuously stimulated its application to cosmology [1]- 
especially to cosmic inflation. Given the fact that the de- 
tection of some stringy effects in accelerators is difficult in 
the foreseeable future, string cosmology presumably pro- 
vides the only way to test the viability of string/M-theory 
concretely. Fortunately the recent measurement of the 
Cosmic Microwave Background (CMB) [2] has brought 
the high-precision cosmological dataset by which infla- 
tionary models can be seriously constrained [3]. 

The development in compactification with flux and 
branes has provided examples of dynamics that fix all 
moduli [4,5]. In particular, the authors of Ref. [5] showed 
that by adding anti-D3 branes to the solutions of [4] , one 
can lift the supersymmetric vacuum energy and find lo- 
cally stable minima with positive cosmological constant 
(KKLT vacua). 

Lately we proposed a string-inspired cosmological 
model [6] based on a massive Dirac-Born-Infeld (DBI) [7] 
scalar field <f> rolling on the anti-D3 brane of KKLT vacua 
(see also Ref. [8]). We showed that this scenario satisfies 
observational constraints coming from the CMB temper- 
ature anisotropics by evaluating the spectra of scalar and 
tensor perturbations generated during the rolling scalar 
inflation [6]. The problem of reheating associated with 
the DBI action [9] is overcome by taking into account a 
negative cosmological constant of the KKLT vacua. The 
origin of dark energy can be explained if the potential 
energy does not exactly cancel with the negative cosmo- 
logical constant at the potential minimum. 

In this work we shall consider several important cos- 
mological aspects which distinguish our scenario from 
others. It is well known that electromagnetic fields are 
not generated in a Friedmann-Robertson- Walker (FRW) 
background if the underlying theory is conformally in- 
variant as in the classical electrodynamics. On the other 



hand we observe magnetic fields with magnitudes of or- 
der 10" 6 G on scales larger than 10 kpc [10]. A number of 
people tried to solve this discrepancy by breaking confor- 
mal invariance of the theory [11] or by breaking confor- 
mal flatness of the background geometry [12,13]. In spite 
of these attempts it is fair to say that it is not easy to 
construct a satisfactory cosmological model based on par- 
ticle physics which explains the origin of seed magnetic 
fields and also satisfies other observational constraints. 
Presumably a unique example known earlier, is provided 
by Gasperini et al. [14] who highlighted the coupling be- 
tween dilaton and gauge fields based upon the Pre-Big- 
Bang scenario. This model respects conformal invariance 
which allows the photon to decouple from the metric. 

The DBI action in our string-inspired scenario involves 
the coupling between gauge fields and the inflaton field 
<j>. We show that this coupling naturally leads to the pro- 
duction of cosmological magnetic fields during inflation 
and the amplitude can be greater than the limiting value 
of the seed fields for galactic dynamo. We also study the 
reheating dynamics in which radiation is generated very 
efficiently through the amplification of gauge fields. 

Model - Let us consider a massive open string exci- 
tation of the anti-D3 brane of KKLT vacua as a candi- 
date of the inflaton <f>. In addition to the massive infla- 
ton potential V((f>), we implement a negative cosmolog- 
ical constant (—A) which comes from the stabilization 
of the modulus fields [4,5]. Then the DBI type effective 
4-dimensional action for our system is described by [6] 



S = J d 4 a;| % /^ (M%R/2 + A) 



-V{4>)^- de%^ + duM^> + F^)\ , (1) 
where V(<p) = f3 2 T 3 e^ m2 ^ , M p is the reduced Planck 



1 



mass, and F„ 



2V[ M A ly ] is the Maxwell tensor with 



the four-potential. Here T% is a brane tension, (3 is a warp 
factor, and m is the mass of scalar vertex operators [6]. 

In a spatially flat FRW metric with a scale factor a, 
we get the background equations of motion 



H 2 = 



1 



3M| 



0/(1- 4> 2 ) + 'iH(t> + (3m 2 



b 2 -A 
6 = 0, 



(2) 
(3) 



where H = a/a is the Hubble rate and a dot denotes the 
derivative in terms of cosmic time t. Here we dropped the 
contribution from gauge fields. During inflation the neg- 
ative cosmological constant is negligible relative to V((f>) 
and inflation is realized for <p 2 < 2/3. The presence of a 
negative cosmological constant leads to a successful re- 
heating in which the energy density of <f) scales as a pres- 
sureless dust. In what follows we consider the generation 
of electromagnetic fields during inflation and reheating. 

Cosmological magnetic fields generated in inflation - 
The action to the second order of gauge field becomes 



A 2 



VA-VA 



1 - , 



(4) 



term in Eq. (6) is approximately given by —f(rj) ~ 
-k 2 <p 2 + {2>/2)(f> 2 (aH) 2 ~ <p 2 /i] 2 under a slow-roll ap- 
proximation. We shall consider a situation in which the 
field <j) rapidly decays to radiation after inflation and the 
gauge field becomes effectively massless during the sub- 
sequent radiation dominant era. Then one can find the 
spectrum of super-Hubble gauge fields by matching two 
solutions between inflation and radiation (see sec. 4.2 of 
[15]). This gives the spectrum V Ak = k 3 / (2ir 2 )\A 2 k \ cx k n 

with n of order (f> 2 [see Eq. (90) of [15]]. Since 4> 2 is much 
smaller than 1 (see Fig. 1), the resulting spectrum is 
nearly scale-invariant. In the case of Z-boson with mass 
M z the spectral index is n = 2(M Z /H) 2 [15]. There ex- 
ists an interesting mechanism to amplify magnetic fields 
by parametric resonance during reheating using a cou- 
pling between a gauge field and a Klein- Gordon scalar 
[16]. In this case the effective mass of gauge fields dur- 
ing inflation is typically required to be larger than H 
for parametric resonance to be efficient, thereby giving a 
blue-spectrum n > 1. 

The energy density stored in a magnetic field mode 
|.B fe | is given by [13] 



PB = \B k \ 2 /(87r) =u 4 \/3 k \ 



(10) 



where we used the coulomb gauge (A = = V • A) , and 
A = A. From this Lagrangian we get the equation for 
the each Fourier component of the gauge field 

A k + H{\ - 30 2 )i fe + (1 - 4> 2 )(k 2 /a 2 )A k = , (5) 

where A; is a comoving wavenumber. Introducing a new 



quantity A k = b x l 2 A k with b = V(<f))/\J 1 - </> 2 , we find 

A' k '+[k 2 -f(r 1 )}A k = 0, (6) 

where a prime denotes the derivative with respect to a 
conformal time r\ = J a _1 di and 

/(„=,V-3«^ + ^ + f--^).(7, 

In an asymptotic past (t] — ► — oo) with </> ^ 0, f(rf) 
is vanishingly small compared to k 2 and the solution is 



given by A\ 



lkr, /V2k. The <f> term gradually be- 



comes important during inflation, which leads to the vari- 
ation of A k . The equation (6) has the following solution: 



A k = a k A\ + (3 k Al* , 
where the Bogolyubov coefficient fi k is 



(8) 



(9) 



The spectrum of the gauge field A k can be evalu- 
ated analytically by using the method in Ref. [15]. Af- 
ter the Hubble radius crossing (k < aH), the mass 



where u> = k/a. In the presence of the f(rj) term, \/3 k \ 2 is 
not generally zero, which corresponds the generation of 
electromagnetic fields. Making use of Eqs. (9) and (10) 
the amplitude of the magnetic field is expressed as 



\B k 



^V l/3fc| 



2tt- 



Vf 



Qif(v)Qkdv 



(11) 



where we introduced a dimcnsionlcss quantity Q k — 
\/2kA k . Note that 77/ is the time after which the particle 
creation ceases. We are interested in the amplitude |-Bf ec | 
on a scale corresponding to the time at decoupling, i.e., 
oJdec = kdec/o-dec — 10 -33 GeV. From the above equation 
we obtain 



\Bt ec \ 



= V27ruj dec H e (a e /a dec )a . 



where a is a dimcnsionless quantity, given by 



11 f (aH) 2 
a e tl e 



(12) 



(13) 



Here g{rf) = (aH)~ 2 f(-n) and the subscript ! e' denotes 
the value at the end of inflation. 

Let us estimate the amplitude of cosmological mag- 
netic fields generated in our scenario. The slow-roll pa- 
rameter e = —H/H 2 becomes greater than 1 when <p 2 
reaches 4> 2 ~ 2/3, from which one can find the value of 
4> at the end of inflation as <f) e = 6.715 /(^m) by us- 
ing Eqs. (2) and (3). The COBE normalization gives 
(3 ~ 10~ 9 for an exponential potential [6], thereby yield- 
ing H e ~ 6.0 x 10" 5 M P ~ 10 14 GeV. 
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FIG. 1. The evolution of a and c/> 2 during inflation for 
f3 — 1CF 9 as a function of the e-folds N = j Hdt. We choose 
the cosmological mode which crossed the Hubble radius about 
56 e-folds before the end of inflation. We find that a reaches 
of order 0.01 at the end of inflation. 

The ratio a e /adec depends on the details of reheating. 
For example if the energy density of the field <f> at the 
end of inflation is converted to almost instantaneously 
to radiation, pR = (it 2 /30)g»T^, where g* <~ 100 is the 
number of relativistic degree of freedom, then we obtain 
the reheating temperature Tr ~ 10 15 GeV for /? = 10~ 9 , 
which yields a e /ad ec ~ T dec /T R <~ 10~ 25 . When the 
thermalization process is delayed, one gets the smaller 
reheating temperature, thereby giving larger values of 
Tdec/TR. By Eq. (12) the amplitude of the cosmological 
magnetic field at decoupling epoch is estimated as 

\B d k ec \ ~ (a e /a dec )a G, (14) 

where we used 1 G 2 /8tt = 1.9089 x 10~ 40 GeV 4 . 

The size of a acquired during inflation is crucially im- 
portant to estimate the amplitude of magnetic fields. As 
shown in Fig. 1, a rapidly grows during the inflationary 
stage. This comes from the fact that the second term 
in Eq. (7) increases due to the presence of the a 2 term, 
which becomes dominant after the Hubble radius cross- 
ing even if the fa 2 term is suppressed during inflation. We 
choose the initial condition fa = in our numerical sim- 
ulations, but fa 2 is of order 0.01 in most stage of inflation 
as seen in Fig. 1 . The evolution of a becomes mild during 
the reheating phase and the particle creation ends when 
the system approaches another asymptotic limit rjf — > oo 
characterized by <ft — > 0. 

In the case of instant reheating with (3 = 10~ 9 , one has 
a e /a dec ~ 10- 25 , yielding \Bf e c \ ~ 10~ 27 G for a = 0.01. 
It is known that seed fields with the amplitude larger 
than |-Bf ec | = 10~ 23 G is required for the galactic dy- 
namo for a flat universe without a cosmological constant, 
but this limit is relaxed up to \B% ec \ = 10~ 30 G in the 
presence of a cosmological constant with tt\ = 0.7 [17]. 
Therefore the value |-Bfc 6C | ~ 10~ 27 G is greater than 
this limiting value. Thus our scenario provides a sound 



mechanism to generate seeds magnetic fields through the 
breaking of conformal invariancc. Larger magnetic fields 
can be obtained if the ratio a e /a de c is higher and this 
depends on the details of reheating. 

Reheating from electromagnetic fields - Let us next 
consider how the gauge fields are produced in reheating 
stage after inflation. We denote the comoving wavenum- 
ber corresponding to the Hubble radius at the end of 
inflation as k — a e H e . Then the energy density of pro- 
duced gauge fields can be expressed as 



Ps{k) = 




One can estimate the size of gauge fields generated dur- 
ing inflation for the mode k <~ a e H e . Since the frequency 
of Ak is approximately given by ui 2 ~ k 2 (l — fa 2 ) for these 
modes during inflation, the Bogolyubov coefficient (3k at 
the end of inflation is estimated as 

/?, ~ (l/4)(0 2 ) (e- 2 ^ - e - 2 ^) , (16) 

where we used Af. ~ A 1 ,. Note that {fa 2 ) is the average 
value of (j> 2 , which means that \/3k\ 2 <C 1 at the end of 
inflation for the mode k <~ a e H e . Therefore ps{k) is 
much smaller than the energy density of the field <fi (p^) 
at the beginning of reheating. 

When the reheating stage starts, the fa 2 term rapidly 
grows to unity as seen in Fig. 1. We numerically found 
that <f> oscillates between — 1 and 1 with a slow adiabatic 
damping due to cosmic expansion. During the oscilla- 
tion, <j) tends to stay around <f> ~ ±1 for some moment of 
time as is checked from Eq. (3). While this oscillation is 
not sinusoidal, one may expect that gauge fields are gen- 
erated by parametric resonance due to the non-adiabatic 
change of the frequency in Ak ■ 

Actually this happens depending on the size of the mo- 
mentum k. The gauge fields Ak for the modes deep in- 
side the Hubble radius (k ^> a e H e ) exhibit rapid growth, 
whereas the modes corresponding to k < a e H e are not 
amplified. This can be understood as follows. Introduc- 
ing a new quantity, Ak — a 1 ! 2 Ak, the equation for gauge 
fields for the mode k » a e H e is approximately given as 

A k + (k 2 /a 2 )(l-^ 2 )Ak~0, (17) 

which has an oscillating term fa 2 . The approximate time 
scale for the oscillation of the field <f> is m~g with m e s = 
\fi3m. We shall define the effective resonance parameter 

k 2 If k \ 2 / a e \ 2 ( H e \ 2 

q = Ja^nJ s = 4 \a~Ji e ) J {^) ' (18) 

which is analogous to the one defined in Ref. [18]. 

The parameter q characterises the strength of paramet- 
ric resonance. The modes deep inside the Hubble radius 
(k 3> a e H e ) correspond to the large resonance parameter 
(?> 1) at the of inflation (here we used the fact that 
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H e is the same order as m e ff). As shown in Fig. 2 we 
find that the production of gauge fields for the modes 
k > 10a e H e is so efficient that its energy density can 
surpass that of the field <p. Note that there is no amplifi- 
cation of the fluctuation 5<fr in our scenario at the linear 
level [6] . The excitation of gauge fields continues until the 
coherent oscillation of the field <fi is broken by the back- 
reaction effect of produced particles or the parameter q 
drops down to of order unity. For the modes k < a e H e 
the resonance parameter q is smaller than unity, which 
means that the particle creation is not efficient as seen 
in Fig. 2. The cosmological magnetic field is involved 
in this case (k <C a e H e ), implying that its generation is 
weak during reheating. 
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FIG. 2. The evolution of ps(fc) during reheating with 
(3 = 1CT 9 for (a) k/(a e H e ) = 1, (b) k/{a e H e ) = 10 and 
(c) k/(a e H e ) = 20. We also show the evolution of the energy 
density for the field <j) (denoted by p^). We find that ps(fc) 
surpasses p$ for the mode k/(a e H e ) > 10. 

In Fig 2 we do not account for backreaction and rescat- 
tcring effects of created particles. This is particularly im- 
portant when ps(fc) catches up p^. In addition the effect 
of conductivity is not negligible if charged particles are 
produced after inflation. This typically suppresses the 
growth of gauge fields [13], whose effect can be impor- 
tant after the intermediate stage of reheating. While the 
explosive production of gauge fields at the initial stage is 
robust, we require a fully nonlinear analysis including the 
effect of conductivity in order to understand the reheat- 
ing dynamics completely. The process of the thermaliza- 
tion is also dependent on the details of such a nonlinear 
stage. It is of interest to do a detailed analysis in such 
a regime, since the reheating temperature is relevant to 
estimate the amplitude of cosmological magnetic fields. 

Conclusions - Generally it is difficult to construct a vi- 
able cosmological model in string theory which satisfies 
all cosmological/observational constraints including in- 
flation, reheating, dark energy and primordial magnetic 
fields. However our scenario using a massive rolling scalar 
on the anti-D3 brane of KKLT vacua can provide a sat- 



isfactory explanation for the above requirements. 

The crucial point for the generation of cosmological 
magnetic fields in our model is that the DBI type action 
breaks conformal invariance due to the coupling between 
gauge fields and the inflaton. This also happens for the 
rolling tachyon field [19], implying that magnetic fields 
can be generated during tachyon inflation. Nevertheless 
tachyon inflation generally suffers from the problems as- 
sociated with large density perturbations and reheating 
[9] . In our scenario these problems are overcome by con- 
sidering a warp metric with a small parameter /3(<C 1) 
and also by a negative cosmological constant, both ap- 
pearing in the KKLT vacua [5] . 

During reheating we find that radiation is efficiently 
produced by the amplification of gauge fields on sub- 
Hubble scales. This corresponds to the "preheating" 
stage as in the bosonic particle creation through the cou- 
pling (l/2)g 2 (f> 2 x 2 [18]. The important point is that we 
used only the DBI type string effective action without 
adding any phcnomcnological terms. It is of interest to 
extend our analysis to the production of fermions using 
the couplings appearing in string theory. 
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